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-<^ i Abstract. In this paper we extend Newton- Steffenssen method for 

jy-^ I solving nonhnear equations, introduced by Sharma [J.R. Sharma, A 

CN ' composite third order Newton- Steffenssen method for solving nonlin- 

ear equations, Appl. Math. Comput. 169 (2005), 242-246] by using the 
-^ I Gejji-Jafari decomposition technique. Several numerical examples are 

^ I given to illustrate the efficiency and performance of this new method. 
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OO ! 1. Introduction 

_il ■ Solving nonlinear equations is one of the most important problems 

O ' in numerical analysis. To solve nonlinear equations, iterative methods 

^^ ' such as Newton's method are usually used. Throughout this paper we 

consider iterative methods to find a simple root a, of a nonlinear equa- 
tion f{x) = 0, where f : I G R ^ R for an open interval I. Many 
/\ . variants of the Newton's method have been suggested in the litera- 

cy ! ture by different techniques. One of them is Adomain decomposition 

method which is used in [B] and other literatures. To implement Ado- 
main decomposition method, one has to calculate Adomain polynomial, 
which is another difficult task. Other techniques have also their limi- 
tations. To overcome these difficulties, a new decomposition technique 
is introduced by Gejji-Jafari in [1]. In this paper we use this technique 
to extend Newton- Steffenssen method introduced by Sharma [2]. 

2. Gejji-Jafari decomposition method 

Consider the nonlinear equation 

fix) = 0. (2.1) 
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Throughout the paper we assume that f{x) has a simple root at a and 
7 is an initial guess close to a. Let us transform the nonlinear equation 
(12.11) into the following canonical form: 

x = c + N{x), (2.2) 

where N{x) nonlinear operator and c is a constant. The main idea of 
this technique is to look for a solution having the series form 



oo 
X 

i=0 



Y.X,. (2.3) 

The nonlinear operator N can be decomposed as 



j~i 



From equations (12. 3p and (12.41) . equation (12.21) is equivalent to 

^x, = c + iV(a;o) + ^<^ivK^xJ - N [Y^xA \. (2.5) 
Thus we have the following recurrence relation: 

Xq = C, 

xi = N{xo), 

X2 = N{xo + xi) - N{xo), 

X2 = N{xo + xi+ X2) - N{xo + xi), 



Xm+l = N{xo + Xi + ... +Xm) - N{xo + Xi + ... + Xm-l)- {'^■^) 

Then 

xo + xi + ... + Xm+l = N{xo + xi + ... + Xm)] m = 1, 2, ..., (2.7) 
and 



00 
X = c 

i=l 

^00 



+ Y.X. (2.8) 



In [T] it is proved that the series Yll^o ^i converges absolutely and 
uniformly to a unique solution of (12. 2p . 
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3. Extension of Newton-Steffenssen method 
Consider the following coupled system: 

9{x) 



fii) + (a; - 7) 



/'(7) + 



0, 



g{x) = f{x)-f{^)-f'{^){x-j). 

The equation (13. ip of the above system can be rewritten as 

/(7)(x-7) 
f'h){x-l)+9{x)' 
Comparing Equations fl2.2p . 1st of fl2.6p and (13.30 . we have 

xo = c = 7. 

and 

/(7)(a;-7) 



f {V{x - 1) + 9{x) 
Note that x is approximated by 

Xm = Xo + 2:1 + ... + Xm, 

where lim^^oo Xm = x. 
For m = 0, 

X^Xq = Xq = C=J. 

For 171 = 1, 

X ^ Xi = xq + xi = 'J + N{xo). 
where A^(a;o) is to be calculate. From (13.51) we have 

f'il)ixo--f)+gixo)' 
Thus (13. 8p becomes 

X ~ Ai — Xo + Xi — 7 - — T — Xo - — ^ r. 1,0. iU 

f ivixo - V + 9{xo) f'{xo) 

which yields the famous Newton's method with second order conver 
gence 

Xn+1 — Xn -777 r. lo.ii 

f [Xn) 

For m = 2, 

X ~ X2 = Xo + Xi + X2 = 7 + A^(xo + Xi). (3.12 
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where A^(xo + xi) is to be calculate. From fl3.5p . (13. 2 p and f IS.lOp we 
have 

N{xo + xi) = ^]/' -^. (3.13) 

Thus we have 

/(7)2 
/ (7)1/(7) - f{xo + Xi)} 

which gives the following well known Newton-Steffenssen method [2] 
with third order convergence 

where Vn = Xn - jr^- 
For m = 3, 

X ^ X^ = xo + xi + X2 + X3 = 'y + N{xo + 0:1 + ^2). (3.16) 

where N{xo + 2:1 + X2) is to be calculate. From (13.51) . (13. 2p and (13.141) 
we have 

/(7)3 
Nixo + X, + X2) = /.(7){/(7) _ /(,„ + ,;)|{/(,„ + ,, + ,,) _ /(7)f -17) 

Thus we have 

X ^ X3 = Xo + Xi + X2 + X3 

^ '^ ~ /'(7){/(7) - f{xo + xSufil) - f{xo + xi + x^^l^^ 
which suggests the following three-step iterative method 



Un •^n 



Xr. 



f{Xnf 
f'{Xn){f{Xn) - fiVn)}' 



"^"^^ " ^'^ f'{Xn){f{Xn)-f{yn)}{f{Xn)-f{Zn)y ^ ^ ^ 
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Similarly we can obtain higher-order iterative methods. For general n 
it can be shown that the (n — l)-step iterative method is 

al = X - ^^"^"^ 






fix.) ' 

fiXn? 



/'(x„){/(x„)-/(al„)}' 



f'{xn){f{xn) - /(al„)}{/(x„) - /(a2„)} 



a(n - 2)„, = x„. - 



f{Xn) 



\n-2 



f'{xn){f{xn) - /(al„)}{/(x„) - /(a2„)}...{/(x„) - f{a{n - 3)„)} 

fixuT-' ^ 

""" f'{xr.){f{xn) - /(al„)}{/(x„) - /(a2„)}...{/(x„) - fiain - hf, 

Now we prove that order of convergence of the iterative method fl3.19p 
is four, which is shown by the following theorem: 

Theorem 3.1. Let a & I be a simple zero of a sufficiently differentiahle 
function f : I C ^ ^ ^ in an open interval I. If Xq is sufficiently close 
to a, then the three-step method defined by fl3.19p has order fourth- order 
convergence. 

Proof. By applying the Taylor series expansion theorem and taking 
account /(a) = 0, we can write 

f{xn) = en + C2el + C3el + C4el + Cr,el + ceel+C7el+csel + 0{el), (3.21) 

where Ck = | ^" , A; = 1, 2, ... and e^ be the error in Xn after n iterations 
I.e. Cfi Xji u^, 

f{xn) = f'{a)[l + 2c2en + 3c,el + Ac^el + 5clei 

+6c6e^ + 7crel + Scge^ + 0(e^)]. (3.22) 

By considering the above relations, one can obtain 

yn = a + C2el + 2(c3 - cl)el + + O(e^). (3.23) 

At this time, we should expand f{yn) around the exact a root by taking 
into consideration (I3.23p . Accordingly, we have 

fiVn) = f\a)[c2el + 2{~cl + c,)el + + 0{e'% (3.24) 
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(3.25) 



From ( Km . (I322D and I^M) it can found that 

Zn = a + cle^ + (3c2c3 - 3c^)e^ + ... + O(e^). 
Now expand f{zn) around the exact a root, we have 

/(z„) = f'{a)[cle' + (3c2c3 - 3c^)e^ + ... + 0(e«)]. (3.26) 

Finally by virtue of ^M>, Km . KM and Km it can be obtained 



that 



e„+i = c^e^ + 0(e^). 



(3.27) 

n 

Similarly we can prove that iterative method (I3.20p has n*'^-order 
convergence. 

4. Numerical Testing 

Here we consider, the following eight test functions to illustrate the 
accuracy of new iterative method. Some of them are taken from [7] and 
some from [8]. The root of each nonlinear test function is also listed. 
All the computations reported here we have done using Mathematica 
8. Scientific computations in many branches of science and technology 
demand very high precision degree of numerical precision. We con- 
sider the number of decimal places as follows: 10000 digits floating 
point (SetAccuraccy=10000) with SetAccuraccy Command. In exam- 
ples considered in this article, the stopping criterion is the |/(a;„)| < e 
where e = 10"^°°*^''. The test non-linear functions are listed in Table-1. 

Here we comparer performance of our new method (I3.19P to the 
methods of Yun (YN) [3], of Chun (CN) |6j and Noor(NR) [5]. The 
results of comparison for the test function are provided in the Table 2. 
It can be seen that the resulting method from our class are accurate 
and efficient in terms of number of accurate decimal places to find the 
roots after some iterations. 

Table 1. Test functions and their roots. 



Non-linear function 




Roots 


fi{x) = sin"^ X — x"^ + 1 




1.4044916482153412260 


f2{x) = 3;2-e^-3a; + 2 




0.25753028543986076046 


hix) = (x - 1)3 - 1 




2 


/4(x) = X^ - 10 




2.1544346900318837218 


hix) = xe^ — sin^ x + 3cosx + 5 


-1.2076478271309189270 


/e(a;) = e^'+7^-30 _ ^ 




3 


//(x) = x^ + sinx + X 







fsix) = sin(2cosx) — 1 — 


^2^gsmx3 


1.3061752018468278250 
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Table 2. Comparison of different methods witli tlie 
same total number of function evaluations (TNFE = 24) 



TestFunction 


Guess 


CH 


YN 


NR 


f09f 


l/ll 


-1.0 


0.15403e-267 


0.41715e-642 


0.52856e-2 


0.22708e-1576 




2.0 


0.35940e-1457 


0.23037e-1530 


0.11251e-3 


0.13526e-2046 




1.0 


0.15403e-267 


0.41715C-642 


0.52856e-2 


0.22708e-1576 


1/2 1 


2.0 


0.74095e-1843 


0.42734C-1881 


0.16354e-4 


0.62927e-2949 




2.5 


0.70566e-1411 


0.31383C-1450 


0.199813-3 


0.18307e-1924 




-1.5 


0.91945e-2055 


0.25223G-2127 


0.46951e-4 


0.13343e-2794 


1/3 1 


3.5 


0.95901e-373 


0.60165e-400 


0.43860e-l 


0.20884e-671 




3.1 


0.22428e-546 


0.32703e-582 


0.13254e-l 


0.36490e-921 




1.5 


0.39254e-Hl 


0.21668e+l 


0.23818e-hl 


0.38950e-663 


I/4I 


1.5 


0.610763-350 


0.25884e-740 


0.27421e-l 


0.47059e-1750 




1.2 


0.38918e-2 


0.12361e-298 


0.85122C+0 


0.716731e-925 




1.0 


0.749233+3 


0.27646e+5 


0.42601C+1 


0.18324e-518 


1/5 1 


-2.0 


0.74075e-140 


0.16833e-150 


0.10709e+l 


0.60969e-322 




-1.5 


0.19698e-1039 


0.21190e-1089 


0.26108e-2 


0.52642e-1607 




-1.0 


0.19824e-693 


0.10164e-1074 


0.38167e-2 


0.49335e-2031 


I/el 


3.5 


0.45434e-5 


0.99135C-6 


0.13063e-|-2 


0.15256e-22 




3.2 


0.11804e-200 


0.14326e-215 


0.92374e-l 


0.69853e-422 




2.9 


0.20627e-l-114 


Indeterminate 


0.42942e-|-0 


0.42509e-496 


IM 


0.3 


0.51508e-2921 


0.46046e-3026 


0.14553e-6 


0.13490e-3702 




0.1 


0.10912e-4558 


0.46696e-4679 


0.77208e-10 


0.16080e-5443 




-0.2 


0.27239e-2693 


0.16105e-2868 


0.19173e-6 


0.17394e-3843 


l/sl 


1.35 


0.29081e-4253 


0.15643e-4384 


0.62249e-9 


0.15434e-5050 




1.31 


0.37148e-8008 


0.11024e-8139 


0.20749e-16 


0.66987e-8948 




1.29 


0.56923e-5177 


0.37152e-5314 


0.56633e-ll 


0.20166e-6195 
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